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In this paper, a detailed analysis of anisotropic effects on the phonon induced spin relaxation rate 
in III-V semiconductor quantum dots (QDs) is carried out. We show that the accidental degeneracy 
due to level crossing between the first and second excited states of opposite electron spin states in 
both isotropic and anisotropic QDs can be manipulated with the application of externally applied 
gate potentials. In particular, anisotropic gate potentials enhance the phonon mediated spin-flip 
rate and reduce the cusp-like structure to lower magnetic fields, in addition to the lower QDs radii 
in III-V semiconductor QDs. In InAs QDs, only the Rashba spin-orbit coupling contributes to the 
phonon induced spin relaxation rate. However, for GaAs QDs, the Rashba spin-orbit coupling has 
a contribution near the accidental degeneracy point and the Dresselhaus spin-orbit coupling has 
a contribution below and above the accidental degeneracy point in the manipulation of phonon 
induced spin relaxation rates in QDs. 
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I. INTRODUCTION 

The study of electron spin states in zero dimensional 
semiconductor nanostructures such as QDs is important 
for the development of next generation electronic devices 
such as spin transistors, spin niters, spin memory devices 
and quantum logic gates.— ~— The electron spin states in 
QDs are brought in resonance or out of resonance by 
applying suitable gate potentials in order to read out the 
spin statesi^ Progress in nanotechnology has made it 
possible to fabricate gated quantum dots with desirable 
optoelectronic and spin properties Very recently, it 
was shown that the electron spin states in gated quantum 
dots can be measured in the presence of magnetic fields 
along arbitrary directions J2r— 

A critical ingredient for the design of robust spin- 
tronic devices is the accurate estimation of the spin re- 
laxation rate. Recent studies by authors in Refs. and 
[B| have measured long spin relaxation times of 0.85 ms 
in GaAs QDs by pulsed relaxation rate measurements 
and 20 ms in InGaAs QDs by optical orientation mea- 
surements. These experimental studies in QDs confirm 
that the manipulation of spin-flip rate by spin-orbit cou- 
pling with respect to the environment is important for 
the design of robust spintronics logic devices J2r— The 
spin-orbit coupling is mainly dominated by the Rashba 20 
and the linear Dresselhaus^ terms in solid state QDs. 
The Rashba spin-orbit coupling arises from structural in- 
version asymmetry along the growth direction and the 
Dresselhaus spin-orbit coupling arises from the bulk in- 
version asymmetry of the crystal lattice £2r2£ Recently, 
electric and magnetic fields tunability of the electron spin 
states in gated III-V semiconductor QDs was manipu- 
lated through Rashba and Dresselhaus spin-orbit cou- 
plings iS^22i 

Anisotropic effects induced in the orbital angular mo- 
mentum in QDs suppresses the Lande g-factor towards 
bulk crysta h 25 ' 29 g-factor can be manipulated through 



strong Rashba spin-orbit coupling in InAs QDs 2 ^ and 
trough strong Dresselhaus spin-orbit coupling in GaAs 
QDs. 25 Large anisotropy effects of the spin-orbit interac- 
tion in self-assembled InAs QDs have been recently stud- 
ied experimentally in Ref. Il2l . In this paper, we study 
the phonon induced spin-flip rate of electron spin states 
in both isotropic and anisotropic QDs. Our studies show 
that the Rashba spin-orbit coupling has an appreciable 
contribution to the spin-flip rate in InAs QDs. How- 
ever, the Rashba spin-orbit coupling has a contribution 
to the spin-flip rate in GaAs QDs near the level crossing 
point and the Dresselhaus spin-orbit coupling elsewhere. 
Anisotropic gate potentials, playing an important role in 
the spin-flip rate, can be used to manipulate the acci- 
dental degeneracy due to level crossing and avoided an- 
ticrossing between the electron spin states |0, 0, — > and 
|0, 1, + >• In this paper, we show that the anisotropic 
gate potentials cause also a quenching effect in the orbital 
angular momentum that enhances the phonon mediated 
spin-flip rate and reduces its cusp-like structure to lower 
magnetic fields, in addition to lower QDs radii. 



The paper is organized as follows: In section II, we 
develop a theoretical model and find an analytical ex- 
pression for the energy spectrum of electron spin states 
for anisotropic QDs. In section III, we find the spin re- 
laxation rate of electron spin states for anisotropic and 
isotropic QDs. In section IV, we show that the cusp- 
like structure in spin-flip rate due to accidental degener- 
acy points in QDs can be manipulated to lower magnetic 
fields, in addition to lower QDs radii, with the applica- 
tion of anisotropic gate potentials. Different mechanisms 
of spin-orbit interactions such as Rashba vs. Dresselhaus 
couplings are also discussed in this section. Finally, in 
section V, we summarize our results. 
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FIG. 1. (Color online) Contributions of Rashba and Dres- 
selhaus spin-orbit couplings on the phonon induced spin-flip 
rate vs magnetic fields in InAs QDs. Inset plots shows the 
energy difference vs. magnetic fields near the level crossing 
point and the p-factor vs magnetic fields. Here we choose 
E = I0 5 V/cm, £ = 20 nm and a = b = 1. The material 
constants for InAs QDs are chosen from Refs. [24] and [30j as 
go = -15, m = 0.0239, j R = 110 A 2 , 7d = 130 eVA 3 , eh u = 



0.54 x 10 5 erg/cm, s; 
and p = 5.6670 g/cm 3 . 



4.2 x 10 s cm/s, s t = 2.35 x 10 s cm/s 



II. THEORETICAL MODEL 

We consider 2D anisotropic III-V semiconductor QDs 
in the presence of a magnetic field along the growth 
direction. The total Hamiltonian of an electron in 
anisotropic QDs including spin-orbit interaction can be 
written a s 18 ' 22 ' 24 



H = H xy 



(1) 



where the Hamiltonian H so is associated with the 
Rashba-Dresselhaus spin-orbit couplings and H xy is the 
Hamiltonian of the electron in anisotropic QDs. H xy can 
be written as 



H 



P 2 



1 



^-2^ + 2 mW> * 



by 2 ) + \g ^B<y z B, (2) 



where P = p + eA is the kinetic momentum operator, 
p = —ih(d x ,d y ,0) is the canonical momentum operator, 

A is the vector potential in the asymmetric gauge, m 
is the effective mass of the electron in the conduction 
band, /is is the Bohr magneton, a — (a x ,a y ,a z ) are the 
Pauli spin matrices, cjq = -77 is the parabolic confining 

potential and £q is the radius of the QDs. 

To find the energy spectrum of Hamiltonian ([2]) , it is 
convenient to introduce the canonical transformation of 
position and momentum operator a o 31 i 32 

x 1 =a 1 / 4 x, P l = ^ J7 P x , (3) 



1 



X2 



^b^y, P<2 = TTTT 



b l/4 P V 



(4) 



Also, the assymetric Gauge potential can be written as 



A, 



-x 2 B 



Vb 



1 A*, 



xiB 



Vb 



(5) 



By substituting Eqs. @1|5} into Eq. ©, we get the 
Hamiltonian in the form: 

h = P 2 +a; 2 +3? (P| + x 2 2 )+p ( Xl P 2 - x 2 Pi)+^g VBVzB, 

(6) 



where h = ^H xy , 



a' P = 



2u) a {b/a) 1 ' i 



Also we use the relation mcjQ7 = 1, where 7 2 = 1 



(U) c /U) ) 



and Lo r 



is the cyclotron frequency. 



(v^+v^) 

The energy spectrum of Hamiltonian © can be found 
as follows. First, we need to find the canonical trans- 
formation U of the four-dimensional phase space, P t = 
(Pi, P2, X\, x 2 ) which diagonalizes the quadratic form of 
the Hamiltonian ([5]). To be more specific, Hamiltonian 
Q without Zeeman spin splitting energy can be written 
as 



/i = P*MP, M = 



1 -p/2 

3 p/2 

Op/2 1 

-p/2 3 



(7) 



where t represents the transpose of a vector. The or- 
thogonal unitary matrix U which exactly diagonalizes 
the matrix M can be written as, 



U 



(*+-*-) 




1 



(8) 



where ps± = 3 — 1 ± d and 



b\4 



1 ± 



" ±= 2 



,(9) 



Va± Vb 



1/2 



(10) 



Vb 




(11) 



In the form of rotated operators P' = UP, the Hamilto- 
nian ^ can be written as 



3S_ + 2) 



Pa 



{Py 



y 



'2 



(12) 



The above Hamiltonian represents the superposition of 
two independent harmonic oscillators. The energy spec- 



trum of H xy can be written as 



%+n _ = (n + + n- + 1) hw++(n + - n~) %u)-+-gQiiB(TzB, 



(13) 
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FIG. 2. (Color online) Spin relaxation rate (1/Ti) vs magnetic fields between the states 1 00+ > and |00 — > in InAs QDs. 
Here we chose to = 20nm (the QD radius). As a reference in Fig. [2jiii) , dashed-dotted line represents the spin-flip rate for 
symmetric QDs with a — b — 3. We see that the anisotropic potential enhances the spin-flip rate and reduces the level crossing 
point to lower magnetic field. 
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FIG. 3. (Color online) Spin relaxation rate (1/Ti) vs QDs 
radius between the states 1 00+ > and |00 — > in InAs QDs. 
Here, we choose B = IT and E — 10 4 V/cm (magnetic and 
electric fields, respectively). We see that the anisotropic po- 
tential reduces the level crossing point to lower QDs radius. 



where n± = a±a± are the number operators. Here, a± 
and <Zj_ are usual annihilation ("lowering") and creation 
("raising") operators. In Eq. Q21 we included Zeeman 
spin splitting energy. 

The Hamiltonian associated with the Rashba and 
linear Dresselhaus spin-orbit couplings can be written 
as^i^i 

Hso = ^ {(T x Py - <JyP x ) + ^ {-<7 X P X + OyPy) , (14) 

where the strengths of the Rashba and Dresselhaus spin- 
orbit couplings are characterized by the parameters cxr 



FIG. 4. (Color online) Spin relaxation rate (1/Ti) vs magnetic 
fields between the states 1 00+ > and |00 — > in GaAs QDs. 
Here we chose, QDs radius, £o ~ 32nm and electric field, 
E — W 4 V/ cm. Again, we see that the level crossing point 
reduces to the lower magnetic fields. For GaAs QDs, we use 
the material constants from Refs. [24] and [3fj as go = —0.44, 
m = 0.067, -jr = 4.4 A 2 , <y D = 26 eVA 3 , ehu = 2.34 x 
10" 5 erg/cm, si = 5.14 x 10 5 cm/s, St = 3.03 x 10 5 cm/s and 
p = 5.3176 g/cm 3 

and cud- They are given by 

f2me \ 2 ^ 3 

a R = j R eE, a D = 0.78 7D f^j E 2 /\ (15) 

where 7^ and are the Rashba and Dresselhaus coeffi- 
cients. 

The energy spectrum of the Hamiltonian associated 
with the Rashba and Dresselhaus spin-orbit couplings 
can be written as 
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FIG. 5. (Color online) Contributions of Rashba and Dres- 
selhaus spin-orbit couplings on the phonon induced spin-flip 
rate as a function of magnetic fields in GaAs QDs. Here we 
choose E — W 5 V/cm, £o = 32nm and a = b = 1. We see that 
only Rashba spin-orbit coupling gives the cusp-like structure 
in the spin-flip rate. 



H so = an (1 + i) 6 1/4 k + (s+ - i) a + + b 1/4: K + (s_ + i) a_ + a 1/4 ?y_ (i - s_) a + + a 1/4 r]- (i + s + ) a_ 

+a D (1 + i) a 1/4 K_ (i - s_) a + + a 1/4 K_ (i + s+) a_ + b l/A rj + (-i + s + ) a + + b l/i rj + (i + s_) a_ + H.c, (16) 



where. 

K± = 

v± = 



1 [1 , eBl 

; r % — Ct ± 1 

2(s+-s_) \£ x h 

1 J 1 .eB£ 
2( S+ - S _) \£^ ±Z ~ 



^ ^ , (18) 



where H.c. represents the Hermitian conjugate, I = 
y/h/mil is the hybrid orbital length and 57 = 
\/wq + uj^/A. It is clear that the spin-orbit Hamiltonian 
and the Zeeman spin splitting energy in both isotropic 
and anisotropic QDs obey a selection rule in which the 
orbital angular momentum can change by one quantum. 

At low electric fields and small QDs radii, we treat 
the Hamiltonian associated with the Rashba and linear 
Dresselhaus spin-orbit couplings as a perturbation. Us- 
ing second order perturbation theory, the energy spec- 
trum of the electron spin states in QDs is given by 



hvo_ 



flUJy 



a 2 n <r_ 



hui x + A 



huiy + A 



.(19) 
,(20) 



where, vu± = lo + ± lu z /2, lj z — A/?i is the Zeeman fre- 
quency, A = goUBB, lu x = lo + + u>-, and u> y = cu + — uj-. 
Also, 



l 



2(s+ 



±—a 2 ± + 2a±/3± T — P±\ ,(21) 

S± S T 



r 



1 I 1 2 fl _ 1 .2 



2( S+ - S _) ^-4-2^^-/3^,(22) 



1/4 J 1 ^ ^ 1 



h f./E + y/b 



,(23) 



•(24) 



III. PHONON INDUCED SPIN RELAXATION 

We now turn to the calculation of the phonon induced 
spin relaxation rate between two lowest energy states in 
QDs. Following Ref. |23, the interaction between electron 
and piezo-phonon can be written a a 18 i 34 



(r,*) = 



2 P Vuj c 



Here, p is the crystal mass density, V is the volume of 
the QDs, b^ a creates an acoustic phonon with wave vec- 
tor q and polarization e a , where a = l,t\,ti are chosen 
as one longitudinal and two transverse modes of the in- 
duced phonon in the dots. Also, v4 qce = qiqk^Pijk&aa 
is the amplitude of the electric field created by phonon 
strain, where q = q/q and e/3yfc = ehi4 for i =^ 
k,i jii 7^ k. The polarization directions of the 
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induced phonon are ei = (sin#cos</>,sin$sin0, cos 6), 
e tl = (cos 9 cos </>, cos 9 sin </>, — sin 9) and e t2 = 
(— sinifi, cos </>, 0). Based on the Fermi Golden Rule, the 
phonon induced spin transition rate in the QDs is given 
bySiM 



where si,st are the longitudinal and transverse acoustic 
phonon velocities in QDs. The matrix element M (qa) 
for the spin-flip rate between the Zeeman sublevels with 
the emission of phonon qa has been calculated pertur- 
bativelyi 34 ' 35 As a result, we have: 



1 



2tt 

T 



d 3 o 

- — -3 V \M(qa) \ 2 5{hs a q- e ,o,- 
( 2n ) a =u 



' £ 0,0,+) : 

(26) 



where, 



^ = c(\M x \ 



(27) 



(eh 14 ) 2 (g/.i B B) 3 
35nh p 



4 1 

3 si 



M x = 



M„ 



{is- + 1) Si (hoj x + A) + {-is- + 1) S 3 {hu) x - A) {-is+ + 1) S 2 (Tk^ + A) + {is+ + 1) S 4 (?iu 



A) 



(7iu; x ) 2 - A 5 



lV4 



A : 



(is+ + 1) Si (7iu; x + A) + (-£s + + 1) S 3 (hu x - A) (£s_ - 1) S 2 (Tiu^ + A) + (-is_ - 1) S 4 (^ - A) 



ft 1 /* 



{tkj x ) 
i 



{Tkjy) 



(28) 



,(29) 



,(30) 



2 («+-«-)' 



[a« {(s+ + 1) j9+ + (1 - a+} + a D {(-s_ - t) a_ + (-1 + is+) /?_}] , (31) 



2(s + 
I 



s-Y 



[a R {(s_ -i)/3 + + {l + i s+ ) a+} + a D {{s+ - i) a_ + (1 + /?_}] , (32) 



"3 



2 (*+-*-)' 
I 



2(s+ 



[or {{s+ -i)/3- + (-1 - is-) a-} + a D {{-s- + i) a+ + (1 + is+) (3+}} , (33) 
- [or {(s_ + i) /3_ + (-1 + is+) a-} + a D {{s+ + i) a+ + (-1 + is_) /?+}] . (34) 



r 



In the above expression, we use c = cil x i + 2ctl x t where 



ga I } I ° ga I 



q-'h 
2puj qa 



and g = 



Also, 



ep.o,- 



-eo,o,+ 



is the Lande g-factor. Also, for longi- 
tudinal phonon modes I x i = I y i = || /if 4 and = 0. 
For transverse phonon modes, J X i 
Ixyt = 0. 

For isotropic QDs (a = 6 = 1, s + 
the spin relaxation rate is given by 



I y t = i§/i?4 and 



1 and S- 



1 _ (e/ii 4 ) 2 (gn B B) 3 



Ti 



357r?i p 



1 4 1 

-+~- 1 i V/ 

s, 3 s 



|Afz 



(35) 

where Mr and Mo are the coefficients of matrix element 
associated to the Rashba and Dresselhaus spin-orbit cou- 
pling in QDs and is given by 



Mr. 



M D 



Qr 

V2hn 



1 



1 



H n +1f) 
1 



i 



i 



A 



1 - 



(36) 



(37) 



Since A = g p,sB is negative for GaAs and InAs QDs, it 
means the degeneracy only appears in the Rashba case 



(see 2nd term of Eq. [36]) and is absent in the Dressel- 
haus case. Similarly, the degeneracy only appears in the 
<?-factor for the Rashba spin-orbit coupling. 29 The degen- 
eracy in the Rashba case induces the level crossing point 
and cusp-like structure in the spin-flip rate in QDs. By 
considering second power of A, the spin relaxation rate 
for isotropic QDs is given by 



1 



(ehuf (g/j B B) 3 ( 1 



1 



35ttTi p 
-2 



4 1 

3s| 



2A 2 

h 4 n 4 



t„2 
\ a R 



\2QJ 



(38) 



It can be seen that the spin-flip rate is highly sensitive 
to the effective g-factor of the electron, Zeeman energy, 
hybrid orbital frequency and cyclotron frequency of the 
QDs. 



IV. RESULTS AND DISCUSSIONS 

In Fig. [TJ we investigate the contributions of the 
Rashba and the Dresselhaus spin-orbit couplings on the 
phonon induced spin relaxation rate as a function of mag- 
netic fields in symmetric InAs QDs. Since the strength 
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of the Dresselhaus spin-orbit coupling is much smaller 
than the Rashba spin-orbit coupling = 3.2 at E = 

10 5 V/cm (see Eq. 1X5]) ). only the Rashba spin-orbit cou- 
pling has a major contribution to the phonon induced 
spin- flip rate. The cusp-like structure is absent (see Fig.Q] 
(dashed line)) and the spin- flip rate (1/Ti) is a monotonic 
function of magnetic field (B) for pure Dresselhaus case 
(ocr = 0). We solve the corresponding eigenvalue prob- 
lem with Hamiltonian (fTJ) by applying the exact diagonal- 
ization procedure and the Finite Element Method , 25 ' 36 
obtaining the energy levels. The inset plots show the en- 
ergy difference vs. magnetic field (Fig.QJi)) and effective 
Lande g- factor vs. magnetic field (Fig. QJii)). It can be 
seen that the level crossing point occurs at B = 3.5T 
which is the exact location of the accidental degeneracy 
point in the spin-flip rate either for pure Rashba case 
(ajj — 0) or mixed cases (both and a n p resent). 
Similar results have been discussed in Refs. [22| and [23| 
and we consider these results as a bench mark for further 
investigation of anisotropic orbital effects on the spin-flip 
rate in QDs. 

Fig. E] explores the anisotropic effects on the spin- 
flip rate vs. magnetic fields for the electric fields E = 
10 4 , 5 x 10 4 , 10 5 Vy cm. It can be seen that the enhance- 
ment in the spin-flip rate occurs with the increase in elec- 
tric fields. The accidental degeneracy point in the spin- 
flip rate is not affected by the electric fields which tells us 
that it is purely an orbital effect and is independent of the 
Rashba-Dresselhaus spin-orbit interaction. In Fig. EJii) , 
the accidental degeneracy point is found at the magnetic 
field B = 3.5T. However, this point increases to the 
larger magnetic field B = 6.2T in Fig. 2(ii). The exten- 
sion in the B-field tunability of the spin-flip rate mainly 
occurs due to an increase in the area of the symmetric 
quantum dots. Note that the area of the quantum dots 
in Fig. Hfii) is 9 times larger than the dots in Fig. [2fi). 
We quantify the influence of the anisotropic effects on the 
spin- flip rate in Fig. EJiii). Here we find that the quench- 
ing in the orbital angular momentu m 26 ' 27 enhances the 
spin-flip rate and reduces the accidental degeneracy point 
to lower magnetic fields (B — 5.85T) compared to the 
symmetric quantum dots {B = 6.2T). As a reference, in 
Fig-HJiii), we also plotted the spin- flip rate vs. magnetic 
fields (shown by dashed-dotted line) for symmetric QDs 
(a = b = 3) at E = 10 4 V/cm and £ = 20 nm. Note that 
the area of the isotropic and anisotropic quantum dots in 
Fig.[2jii) and Fig.^iii) are held constant. The expression 
for the level crossing point is given by the conditio n 22 ' 23 

£o.o- = £o,i,+ i- e -> ft ( UJ + - w -) = \9o\vbB (see Eq. [13]). 
For isotropic QDs (a = b = 1), the condition for the level 
crossing point is Q — cu c /2 = \gofisB /%. It means, when 
the difference between the hybrid orbital frequency to the 
half of the cyclotron frequency becomes equal to the Zee- 
man frequency then the degeneracy appears in the energy 
spectrum which give the level crossing point and cusp- 
like structure near the degeneracy in the spin-flip rate in 
QDs. If we compare the condition of the level crossing 
point for the isotropic and anisotropic QDs, we find that 



the anisotropic potential reduces the level crossing point 
to lower magnetic fields if the area of the symmetric and 
asymmetric quantum dots is held constant. However, if 
we increase the area, the level crossing point extends to 
the larger magnetic field. 

In Fig. [3l we study anisotropic effects on the phonon 
induced spin-flip rate vs. QDs radii in InAs QDs. Similar 
to Fig.^iii), the anisotropic potential enhances the spin- 
flip rate and reduces the accidental degeneracy point to 
lower QDs radii at £q = 43nm. 

Next, we investigate the phonon induced spin relax- 
ation rate in GaAs QDs. In Fig. 01 we plot the phonon in- 
duced spin-flip rate vs. magnetic fields for both isotropic 
(a = b = 8, (solid line)) and anisotropic (a = 1,6 = 64 
(dashed line)) GaAs QDs. It can be seen that the cusp- 
like structure due to the accidental degeneracy can be 
manipulated to lower magnetic fields in the phonon in- 
duced spin-flip rate with the application of anisotropic 
gate potentials. 

The contributions of Rashba and Dresselhaus spin- 
orbit couplings on the phonon induced spin-flip rate 
vs. magnetic fields in GaAs QDs are shown in Fig. [5] 
The Rashba and Dresselhaus spin-orbit couplings become 
equal at very large electric field E — 3.02 x 10 6 V/cm in 
GaAs QDs. Below this value of electric field, only the 
Dresselhaus spin-orbit coupling has a major contribution 
on the phonon induced spin-flip rate in GaAs QDs. How- 
ever, near the level crossing point (for example, B = 6T 
in Fig. [SJ , the accidental degeneracy appears due to only 
the Rashba spin-orbit coupling, which gives the cusp-like 
structure in the spin-flip rate. 



V. CONCLUSIONS 

In summary, we have analyzed in detail anisotropy ef- 
fects on the electron spin relaxation rate in InAs and 
GaAs QDs, using realistic parameters. In Fig.[U we have 
shown that only the Rashba spin-orbit coupling has a 
major contribution on the phonon induced spin-flip rate 
in InAs QDs. In Fig. El El and E] we have shown that 
a cusp-like structure due to the accidental degeneracy 
point appears in the phonon induced spin-flip rate and 
can be manipulated to lower magnetic fields, in addition 
to lower QDs radii, with the application of anisotropic 
gate potentials in III-V semiconductor QDs. Also, we 
have shown that the anisotropic gate potential causes a 
quenching effect in the orbital angular momentum that 
enhances the phonon induced spin-flip rate. Finally, in 
Fig. [5] for GaAs quantum dots, we have shown that the 
Dresselhaus spin-orbit coupling has a major contribution 
on the spin-flip rate before and after the accidental de- 
generacy point, and the Rashba spin-orbit coupling has a 
contribution near the cusp-like structure. These studies 
provide important information for the design and control 
of electron spin states in QDs for the purposes of build- 
ing robust electronic devices and developing solid state 
quantum computers. 
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